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We investigate ionization of neon atoms by an isolated attosecond pump pulse in the presence of 
two coherent extreme ultraviolet or x-ray probe fields. The probe fields are tuned to a core-valence 
transition in the residual ion and induce spectral shearing of the photoelectron distributions. We 
show that the photoelectron-ion coincidence signal contains an interference pattern that depends 
on the temporal structure of the attosecond pump pulse and the stimulated core-valence transition. 
Many-body perturbation theory is used to compute “atomic response times” for the processes and 
we find strikingly different behavior for stimulation to the outer-core hole (2 p «-► 2s) and stimulation 
to the inner-core hole (2 p <-► Is). The response time of the inner-core transition is found to be 
comparable to that of state-of-the-art laser-based characterization techniques for attosecond pulses. 


I. INTRODUCTION 

Atoms and molecules are today routinely probed and 
controlled on the atomic time scale in various branches of 
attophysics p]. Tailored laser fields are used to control 
electron trajectories and to probe high-order harmonic 
generation (HHG) [2j [3J. The combination of phase- 
locked attosecond (as) extreme ultraviolet (XUV) pulses 
and femtosecond (fs) infrared (IR) laser fields has found 
numerous applications: The IR-field can act as an intense 
probe to break chemical bonds [4] or map time into fre¬ 
quency space by the so-called attosecond streaking tech¬ 
niques [Mg. Alternatively, the XUV pulse can act as a 
probe to study electron/hole dynamics by transient ab¬ 
sorption techniques mm- XUV and x-ray free-electron 
lasers (FELs) p"2liT5j are accelerator based sources that 
provide pulses of fs duration with peak powers reaching 
the gigawatt (GW) range. The capabilities of these new 
sources are rapidly evolving, including the demonstration 
of wavelength-tunable pulse pairs (T6HT8] and spectral- 
temporal pulse shaping by seeded FELs m, with pos¬ 
sible production of GW-as pulses [20H23] . Diverse fs 
pump-probe schemes at FELs led to groundbreaking ex¬ 
periments in chemical reaction dynamics (24JI27] and the 
extension of these techniques to the as timescale could 
open new avenues for the observation and control of elec¬ 
tron dynamics. Attosecond pulses have mainly been 
characterized using streak-camera techniques EH2, or 
interferometric techniques [28ti30| . where photoelectrons 
are treated as “replicas” of the corresponding attosecond 
pulses that are controlled by a phase-locked optical laser 
field. The implementation of such techniques at FELs 
is cumbersome due to challenges of synchronization of 
FEL pulses to optical lasers (SH E3 Furthermore, re¬ 
cent difference measurements of photoelectrons resulting 
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from different atomic initial states have evidenced that 
such XUV-IR schemes, which rely on laser-driven photo¬ 
electron dynamics, suffer from uncertainties on the order 
of tens of attoseconds 0(33]. 

This article explores the possibilities to exert control of 
a photoelectron without laser-electron continuum inter¬ 
action. To this end, we study a process where an isolated 
XUV-as pulse ejects an electron from neon, while a pair 
of XUV (or x-ray) probe fields drives a transition in the 
ion. The probe fields are tuned to predominately stimu¬ 
late hole transitions rather than electron continuum dy¬ 
namics. A small red and blue detuning of the probe fields 
relative to the resonant hole transition is used to induce 
spectral shearing of the photoelectron replica to lower 
and higher kinetic energies, respectively. The shearing 
process is illustrated under pathway (S+) in Fig. [l] (a) 
relative to (1), the unshifted one-photon ionization pro¬ 
cess. This opens up for a novel control scheme of the fi¬ 
nal electron energy by indirect interaction with the probe 
fields via the hole in the remaining ion. Interestingly, we 
will show that correlation between the electron and the 
hole is not required to explain the mechanism and that 
the process can be described within the independent par¬ 
ticle approximation. 

The paper is organized as follows. In Sec. [II] we in¬ 
troduce the proposed scheme and outline our primary 
theoretical method for the study, which is based on a 
one-dimensional model of the neon atom. In Sec. m 
we present our main res ults fo r the time-dependent one¬ 
dimensional model [Sec. Ill A| and give an interpretation 
of the result using three-dimensional perturbation theory 
[Sec. IIIB . In Sec. |IV| we discuss our results and give an 
outlook for experimental realization of the scheme. Fi¬ 
nally, in Sec. [V] we present our conclusions. 


II. METHOD 

As illustrated in Fig. [I] (a), we consider photoioniza¬ 
tion by an isolated XUV-as pulse in the presence of two 
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FIG. 1, (Color online) (a) Photoionization processes in neon: 
(1) photoelectron replica of attosecond pulse by one-photon 
absorption from 2 p state; and (S+) non-sequential two-photon 
processes generating spectrally sheared replicas by stimulated 
ionic transitions with a final hole in the 2s state, (b) Sketch of 
the proposed experiment where an IR laser field is split into 
two parts for generation of the XUV pump (1) and XUV probe 
(2,3) fields by HHG. Photoionization of neon atoms is then 
studied with electron-ion coincidence detection, (c) Photo- 
electron spectrum with residual hole in 2 p state; (d) photo- 
electron spectrum with residual hole in 2s state, both com¬ 
puted by 1D-TDCIS (within an independent-particle model). 
See main text for the labeled spectroscopic structures. 


coherent XUV (or x-ray) probe fields. The complex XUV 
amplitude of the incoming field on the neon atom written 
as (atomic units, h - e - m - 47reo = 1 , are used unless 
otherwise stated) 


E(cj) — Ei (cj) + F/2 $(c<J — ^2 ) + EqS^CjJ — CJ3 ), (1) 

where the pump amplitude is Ei(uj) 
\Ei(cj)\exp[i(j)i(uj)] with a central frequency uj 1 
and a spectral bandwidth Ac o\. The peak intensity 
of the pump pulse is set to I\ - 7 x 10 12 W/cm 2 with 
a Fourier-limited pulse duration of 244 as. The group 
delay 


Ti(w) 


#1 
dui ’ 


( 2 ) 


describes the arrival time of a particular frequency com¬ 
ponent uj of the pump pulse at the target. We consider 


probe frequencies ojf= 2 , 3 , that are symmetrically red and 
blue shifted relative to the 2 p ++ 2 s hole transition 


<^2 “ e 2p ~ £2 s ~ 

^3 “ ^2 p ~ €2 s + Sid, 


where the outer hole energy is € 2 P = -21.56 eV and the 
inner hole energy is 62 s = -48.47 eV [34]. The peak inten¬ 
sity of the probe fields is set to If = 3.5 x 10 12 W/cm 2 . 
The detuning is supposed small compared to the band¬ 
width of the pump field, Su « Acji, and the probe fields 
are quasi-monochromatic with a bandwidth much smaller 
than the detuning, Aco> 2,3 « Scj, as indicated by the delta 
functions in Eq. |lj) . The group delay of the probe fields 
is defined as 732 = (03 - 02 )/ 2 &j using the spectral phase 
of the probe fields, 02,3 = argjE^s}- In Fig. [l] (b) we 
propose an experimental setup where the pump (1) is 
generated by HHG from an ultra-short IR pulse and the 
probe fields (2,3) are odd HHG harmonics from a longer 
IR pulse. The group delay of the probe fields is then 
locked to the phase of the IR laser field due to the non¬ 
linear HHG process, 732 °c 0 ir, and can be accurately 
controlled by an IR laser-delay stage. Implementation of 
this probe technique at FEL sources would require two- 
color pulse pairs [MU] and accurate temporal and phase 
control [T9] . 

As depicted in Fig jl] (a), the proposed probe process 
can be cast in terms of single-particle transitions. An ap¬ 
propriate numerical method is therefore time-dependent 
configuration interaction singles (TDCIS) [35]. For the 
calculation of the photoelectron spectrum we imple¬ 
mented the coupled surface flux method [36] . similar to 
the approach used in Ref. m To discuss the basic pro¬ 
cess, we opt for a one-dimensional (ID) description of 
the considered process. In Fig. l] (c) and (d) we show 
simulated ionic channel resolved photoelectron spectra 
for a ID model of the neon atom. Surprisingly, we have 
found that it was important to account for stimulated 
ion dynamics by the probe fields after the electron has 
escaped the inner region, as explained in Appendix A. 
Electron-electron correlation effects do not influence the 
ionization process significantly and this allows us to fur¬ 
ther approximate the system by an independent-particle 
model, where the correlation coupling terms in TDCIS 
are neglected. This makes the interpretation of the nu¬ 
merical results more tractable but also helps to speed-up 
the computational time. 


III. RESULTS 


In Sec. Em we present our numerical results for the 
time-dependent ID model of neon. In order to inter¬ 
pret our results and to make quantitative estimates we 
then turn to perturbation theory in Sec. Ell where we 
first co nsider the stimulated 2 p ++ 2 s (XUV) transition 
in Sec. IIIB 1 and then the 2 p ++ Is (x-ray) transition in 
Sec. IIIIB 21 
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Energy of photoelectron (eV) 


FIG. 2. (a)-(c) (Color online) Normalized photoelectron 

distribution of the (S+) peak in Fig. [l] (d), as a function of 
phase difference between probe fields, 032- Data is computed 
using 1D-TDCIS (within an independent-particle model), (a) 
Fourier limited pump pulse; (b) linear chirp; and (c) quadratic 
chirp of pump pulse, (d)-(f) same as (a)-(c) but with nor¬ 
malization at each individual kinetic energy. The left vertical 
axis labels the relative phase in radians, while the right axis 
labels the extracted group delay of the attosecond pulse in 
femtoseconds defined in Eq. shown by the dashed white 
curve. 


FIG. 3. (Color online) Detailed study of the “response time” 
of the (S+) peak for three different detuning of the probe field 
given an unchirped pump pulse (a = /3 = 0). In accordance 
with Eq. 0, the response time T p b, is extracted by making 
a cosine fit to the phase-dependent oscillations of the photo- 
electron probability, e.g. the modulations shown in Fig. [1(d) 
(where Scj = leV). The raw data (not shown) has been fitted 
to a line in order to extract the o-parameter of the pump 
pulse. Standard deviation of the linear fit is indicated by the 
error bars. The data was computed by the ID time-dependent 
independent-particle model discussed in Appendix A. 


A. Time-dependent ID model 


Photoelectrons leaving the residual ion with a hole in 
the 2 p state [Fig. [l] (c)] exhibit one broad peak (1) due 
to absorption of one pump photon (with u)\ = 68 eV and 
Acoq = 7.5 eV) and two narrow peaks (2,3) due to ab¬ 
sorption of either probe field (with ^2,3 = 26.9 =f 1 eV and 
Aco> 2,3 = 0.125eV). The broad peak (S-) is due to ab¬ 
sorption of one pump photon and stimulated emission of 
one probe photon. Weaker peaks labeled with (P 2 ) and 
(P 3 ) denote two and three probe-photon processes, re¬ 
spectively. Photoelectrons leaving the residual ion with 
a hole in the 2 s-state [Fig. IT] (d)] exhibit a peak (P) from 
the pump field and a peak (S+) due to absorption of one 
pump photon and one probe photon. 

In Fig. [2] we show the behavior of the (S+) structure 
as a function of relative probe phase 0 32 = 03 - 02 , for (a) 
a Fourier limited pump pulse, 0 i(cj) =0; (b) a quadratic 
phase dependence, 0 i(cj) = a(uj - uj\) 2 with a = 100; 
and (c) a cubic phase dependence, 0 i(cj) = /3(uj - coq) 3 
with /3 = 100. In Fig. [2] (a) the (S+) peak vanishes 
for 032 « 0 , while the peak is maximized for 0 32 « 7 r. 
As will be derived in Sec. m Bj this is due to a rela¬ 
tive 7r-shift between the two-photon paths that have pos¬ 
itive and negative detuning with respect to the hole res¬ 
onance, respectively. In Fig. [2] (d)-(f) we show more 
clearly the 032 -dependence of the (S+) structure by di¬ 


viding every phase-dependent curve, at a fixed energy of 
Fig. 2. (a)-(c), by its maximal value. In analogy with 
spectral shearing interferometry, the modulations (indi¬ 
cated by white dashed curves) are expected to depend on 
the group-delays of the attosecond pump pulse with (d) 
constant value (e) linear chirp and (f) quadratic chirp, re¬ 
spectively. By direct fit to the unchirped case (a = ft = 0) 
we obtain an extracted a-value of -0.288, corresponding 
to a drift of -0.5as/eV. Where does this “response time” 
come from? 

Further, our simulations show that the (S-) peak 
[Fig. 0 (c)] has a similar 032-dependence as (S+). The 
one-photon peaks ( 1 ) and ( 1 ’) also modulate with 0 32 , 
but the variation is opposite to that of (S+) and (S-). 
Physically, this is because the probe fields are effectively 
shifting the ionic channel of the photoelectrons, e.g. by 
redistribution of population from peak ( 1 ) to peak (S+). 
The total photoelectron spectrum, unresolved for the 
residual ionic state, does not show clear 023 -dependence. 
Unfortunately, this makes the experimental measurement 
of the (S+) modulations challenging, because it must rely 
on coincidence detection, as we will discuss in Sec. m 

In Fig. [3] we show the “response time” for the case of 
an unchirped pump pulse, i.e. by zooming in on the 
dashed curve in Fig. [ 2 ] (d), for three different symmetric 
detuning of the probe fields, = 1 , 1.5 and 2 eV. All de¬ 
tuning show qualitatively the same result with a response 
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a 

6w (eV) 

Awi (eV) 

A u f (eV) 

a 

Numerical 

with ATI 

Without ATI 

Difference 

-100 

1 

7.5 

0.0625 

-99.79 

-100.127 

-99.85 

-0.277 

100 

1 

7.5 

0.0625 

99.79 

99.55 

99.86 

-0.31 

100 

1 

7.5 

0.125 

99.5 

99.058 

99.39 

-0.332 

10 

1 

7.5 

0.125 

9.953 

9.646 

9.943 

-0.297 

100 

0.5 

7.5 

0.125 

97.69 

97.149 

97.57 

-0.421 


TABLE I. Comparison of the retrieved a from numerical calculations and analytical estimate (a). 


time in the range 40-55 attoseconds. All curves exhibit a 
negative slope with extracted a -values for detuning 1 eV, 
1.5 eV and 2eV equal to -0.288, -0.239 and -0.160, re¬ 
spectively. In order to interpret this behavior we have 
additionally performed simulations where the photoelec¬ 
tron does not interact with the probe field. In Fig. [3] 
we label this result as “no ATI” (no above-threshold ion¬ 
ization) because the electron can not absorb probe pho¬ 
tons after it has been ejected within this model (see also 
Appendix A). Interestingly, the extracted (^-parameter is 
-0.00294, which is much closer to the expected zero value. 
This shows that XUV driven electron continuum dynam¬ 
ics must be responsible for the finite response time. The 
discrepancy in recovering the parameters of the pump 
pulse is attributed to an “atomic response time” and a 
derivation of how the group delay ti(cj), is mapped to 
the (S+) structure will be given by perturbation theory 
in the next section. In the case of a chirped pump pulse, 
e.g. a * 0 case, the finite duration of the probe fields will 
also affect the extracted values. In the simplified case 
of Gaussian pulses (with linear chirp given by a) and by 
considering the “no ATI case”, we have found that this 
effect can be approximated as 

/ 2 A UJ 2 n A UJ 2 f \ 

a = a (i | 4l n (2)(5a; 2 ) ’ (3) 

where a is the extracted value for the ^-parameter of 
the pulse. The bandwidths of the pump and probe 
pulses are labeled as Auq and Acjj = 2 ; 3 , respectively, and 
defined as the full-width at half-maximum (FWHM), 
Ei(uj) ~ exp[-21n(2)(cj - c^) 2 /Acjf]. In table [i] we show 
the reasonable agreement between this simple analyti¬ 
cal estimate (a) and the numerical simulations (without 
ATI). Further, we note that the difference between the 
case with ATI and without ATI (right-most column) with 
Slu = leV is approximately -0.3 in agreement our finding 
for the a = 0 case. This indicates that the electron con¬ 
tribution does not depend strongly on the chirp of the 
pump pulse. Finally, the difference between the numeri¬ 
cal simulations with a = 100 and Sco = 0.5 eV (in table [T]) 
shows a larger negative slope than the Scu = leV case, 
in agreement with the trend found for the a = 0 case 
(in Fig. [3]). In the following we will not focus on these 
detailed pulse convolution effects, that occur due to the 
finite bandwidth of the probe fields A cjf, but rather ex¬ 
plain the fundamental reason for why electron continuum 


dynamics leads to a non-zero response time. 

B. Time-independent perturbation theory 

In order to better understand the physical mechanism 
of the atomic response time we now turn to perturbation 
theory. The dominant complex amplitudes that give rise 
to 032 -dependent modulations of the (S+) structure are 
given by 

S p b,f = -rE 1 (id pb -u)f)E f M pb j, f = 2,3 (4) 

corresponding to absorption of one photon from the 
pump field with energy uj[ - uj pb - (Jf, followed by one 
from either probe field with energy ujf. 

1 . Stimulation of 2 p -> 2s hole transition 

In Eq. Q, the two-photon matrix element M pb j = 
M p b(u)pb-Wf,Wf), describes a transition from the ground 
state to a final state with one electron in the continuum 
state p = fcs, kd (with energy e p > 0) and a hole in the 
atomic orbital b - 2s (e b < 0). Energy conservation is 
imposed as uj pb - e p -e b = The probability density 

for electrons in the (S+) structure is computed as the 
square the two complex amplitudes with f - 2,3 leading 
to an interference pattern over 032 °c T 32 

W p b » \S p b,2 + S p b, 3 I 2 

=\A pb \ - | B pb \ cos [26u (n - 732 + r pb )] , (5) 

where \A pb \ is the incoherent sum of the transition 
strengths, while | B pb \ relates to the cross term of the 
amplitudes. We note that the group delays of the pump 
r 1 and the probe fields 732 enter with opposite signs in 
Eq. which must be the case because if the pump 
field is delayed by a certain amount the probe field must 
also be delayed by the same amount to recover same out¬ 
come. Besides the group-velocity delays t\ and 732 , the 
interference pattern is delayed by 

T p b = [arg (M pb:2 M* b 3 ) - 7 r]/ 2 < 5 w, ( 6 ) 

which depends on the phase difference between the two- 
photon matrix elements, and can be interpreted as an 
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atomic response time for creating the (S+) peak. Be¬ 
cause it is convenient to define the response time as a 
small number we have included a -it inside the square 
bracket in Eq. © to remove the relative 7r-shift between 
the two-photon transitions due to the resonance. It fol¬ 
lows from this definition that there is a minus sign on the 
cosine in Eq. (J5|. While may be regarded as a lim¬ 
iting factor for determining the unknown ti, explaining 
the modest error in the extracted <a-value above for the 
ID model, the response time is an interesting quantity to 
study further, as it contains information about the stim¬ 
ulated core-valence transition, in particular, it contains 
information about the phase difference of the two-photon 
(XUV-XUV) matrix elements. 

In order to make a quantitative estimate of we 
now turn to many-body perturbation theory to describe 
neon in 3D and include correlation effects. Following 
Ref. [38] . our calculations are based on single-particle 
states that are expanded on a spherical basis 0*(r) = 
The radial wavefunctions Ri(r), are 
eigenstates to the restricted Hartree-Fock (HF) equation 
for occupied states, while the unoccupied (virtual) states 
are additionally attracted by an effective spherical po¬ 
tential to model the long-range Coulomb interaction be¬ 
tween electron and ion. The two-photon matrix ele¬ 
ments are separated into two terms using second quanti¬ 
zation 


Mn 


pbj 1 


M P bj +M P bj > 


(7) 


because the probe field can either stimulate a hole tran¬ 
sition or a continuum electron transition, as illustrated 
in Fig. [5] (a) and (b), respectively. The stimulated hole 
term 


Tiy-(hole) _ Za'bZpa' 

pbJ 1/ (u f -e a r + e b )’ 


(8) 


describes a dipole transition of an electron from any oc¬ 
cupied single-particle state </v, to the final electron state 
followed by a dipole transition of the hole to the fi¬ 
nal state q V The radial orbitals are chosen to be real, 
which implies a real transition to a given final partial 
wave state, (\> p = ks, kd for = 2s. In contrast, the stim¬ 
ulated electron term is complex 


0 = lim f - --- 

pbJ €- 0 + 2? (w[ - ej/+ e b + 1 £) 


=p.v. 


Zpp' Zp'b 


(wj -e p > +e b ) 


■ 11TZp r Z r ij) 


(9) 


and it describes an initial dipole interaction that excites 
an electron from the occupied state = 2s, to the un¬ 
occupied states, </y = n'p and k'p. The second dipole 
interaction then stimulates an electron transition in the 
continuum +e^ > 0 as shown in Fig. [ 5 ] (c)) to the final 
state 4> p = ks,kd. In Eq. (|9| we write the matrix ele¬ 
ment as real non-resonant contributions (principal-value 
sum-integrals over p' = k'p ) and an imaginary resonant 
contribution (via the intermediate state <fi r = k r p with 



70 80 90 100 110 120 


Photon energy of pump field, [eV] 

FIG. 4. (a) (Color online) Squared matrix elements \M p bj\ 2 , 

for final s-wave (dotted) and d-wave (dashed), including both 
stimulated hole and electron paths in bold line and only elec¬ 
tron path in thin line, (b) Response time for photoemission 
along polarization axis, s-wave and d-wave. Data presented in 
(a) and (b) are computed by a 3D independent-particle model 
of neon, (c) Response time for photoemission along the polar¬ 
ization axis (within a correlated model including both time- 
orderings) for the stimulated valence hole (2 p 2s) and core 
hole (2 p -> Is) transitions. The streak-camera delay from 
initial 2p (2s) state [38] is shown for reference. 


e r = eb + bo[). If the probe field is in the IR range the 
stimulated electron transition is a good approximation 
for the total two-photon matrix element M p bj » M^ c '\ 
but this is not the case for the processes studied here with 
nearly resonant XUV transitions. In Fig. [4] (a) we show 
that stimulated electron contributions \M^ le ^ c ^\ 2 are 2- 
3 orders of magnitude smaller than the total contribu¬ 
tions \M p bj \ 2 that are dominated by the strong reso¬ 
nant coupling in the residual ion. Here, the total ma¬ 
trix element is better approximated as a hole transi¬ 
tion plus a small imaginary electron transition M p bj » 

If we note that = -M^ e) 
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FIG. 5. (a) Two-photon diagram for stimulated hole 

transition, 2 p -> 2s, after photoemission from outer state, 
2 p -> ks , kd. (b) Two-photon diagram for stimulated elec¬ 
tron continuum transition from an inner valence state, 2s -» 
k'p ks, kd. (c) Two-photon diagram for stimulated virtual 
electron transition from a core state, Is -> n p -> ks , kd. 

and assume that then 


electrons with momentum k = kz along the polarization 
axis, computed by the complex final state 

Mr) - E i L e- ir >^Y^ M ^)Rk,L(r)yL,M(r), (11) 

L,M 

where k) are scattering phases of the real radial func¬ 
tions, Rk,L,(r). The angle-resolved emissions has a linear 
drift of 0.634 as/eV, quite close to the dominant d-wave. 
Over a large energy range, from 65 to 120 eV, the de¬ 
viation from this linear fit is less than one attosecond. 
The response times of the 3D calculation are in qualita¬ 
tive agreement with those of the ID case, with a delay 
on the order of tens of attoseconds. Surprisingly, we find 
that the slopes of the response times are different in the 
ID and 3D case. It remains an open question if this dis¬ 
crepancy is entirely due to pulse convolution effects or 
if the different electronic structure between ID and 3D 
plays a role. In order to answer this question it would 
be beneficial to perform 3D TDCIS calculations m, but 
this remains an endeavor beyond the scope of the present 
paper. 

Next, we add correlation effects by implementing the 
random phase approximation with exchange (RPAE) on 
the first dipole interaction [351, which increases the lin¬ 
ear slope to 0.733as/eV (not shown). Including also the 
non-resonant, reversed time-order processes (TO), where 
the probe photon is absorbed before the pump photon 
changes the slope marginally to 0.724 as/eV. Although 
the contribution from the reversed-TO is rather small, 
we stress that the ^-dependence reported for the ID 
model in Fig. [3j can not be explained without taking 
this effect into account. 


Tpb 


arg [M pM ] 


Suj 


= — arctan 
SuJ 


7T ZprZrb 
_ Z'abZpal 


7T Zp r Z r }) 

5 

ZabZpa 


( 10 ) 


where zij > 0 are dipole matrix elements between be¬ 
tween the real single particle states fa and (j)j. The 
last step in Eq. ( |Io| ) is valid for small detuning, Suj « 
^abZpalnZprZrbV Interestingly, Eq. (10) shows that T p b 


does not depend strongly on Suj but rather gives direct 
information about the ratio between dipole matrix ele¬ 
ments of the stimulated electron and hole transitions. In 
other words, decreasing the detuning to enhance further 
stimulate the hole transitions will not alter the response 
time. This has been verified by many-body perturbation 
theory for Suj = 1 and 1.5 eV, where we found that the 
response time changed by less than one attosecond. 

In Fig. |4] (b) we present r p 5 from Eq. (| 6 ) using the 3D 
independent particle model given by Eq. (71). A positive 
linear drift is found on both final partial waves, ks and 
kd , which we attribute to an increasing relative contribu¬ 
tion from the resonant electron path, i.e. the numerator 
on the right side of Eq. |To| . We also show T p b for photo- 


Finally, we note that the ratio of the two-photon 
(pump+probe) and one-photon (pump) transition rates 


R 


RfZba 

Si J 


( 12 ) 


scales inversely with the squared detuning of the probe 
field. This implies a boost of the two-photon transi¬ 
tion by tuning the probe fields closer to the resonance. 
Assuming Suj - leV = 0.0358 au, z a \) « 1 au and 
R = 1%, we estimate the required probe field intensity 
is If\R=i% = 7.4 x 1CT 5 au= 2.6 x 10 12 W/cm 2 . In the case 
studied here with one short pump pulse and two long 
probe fields it is only the temporal overlap of the two 
pulses that will contribute to the two-photon transition. 
Using the time-dependent model, discussed in Sec. m3 
we have verified that the probability for the two-photon 
transition does not depend on the duration of the probe 
fields, but rather on the instantaneous intensity and de¬ 
tuning of the probe fields. We refer the interested reader 
to Ref. |39] for an insightful discussion about pump-probe 
schemes on the attosecond time scale. 



























7 


2. Stimulation of 2p -► Is hole transition 


The squared two-photon matrix element for the 2 p -> 
Is stimulated hole transition is roughly two order of 
magnitude smaller than that of the 2 p -> 2 s, shown in 
Fig. El (a), but the trend is otherwise similar. This is 
easy to understand because the dipole coupling from the 
2 p valence state to the Is inner core state is smaller than 
that of 2 p to 2s. In Fig. [4] (c) we compare r^, includ¬ 
ing correlation and both time-orders (RPAE+TO), for 
XUV-stimulated outer-core-valence transition ( 2 p -> 2s) 
and x-ray inner-core-valence transition ( 2 p -> Is). Evi¬ 
dently, the core transition has a much shorter response 
time. This can be explained by the fact that the (main) 
electron path no longer goes through the continuum, but 
instead on a virtual bound excitation 


M 


(elec.) 


<i' Z n 'b 


pb,f 


V (^1 ~ e n' +e&)’ 


(13) 


where uj[ + < 0 for b = Is, as illustrated in Fig. [5] (c). 

In the present calculation we used the real HF energy 
for the Is orbital = -891.70 eV, which should be 

valid provided that the pump pulse is sufficiently short. 
Nonetheless, we have also tested to give the ls-energy an 
imaginary part (equal to 0.27 eV) to mimic the decay of 
the core hole, but this did not change the response time 
by more than one attosecond. 


IV. DISCUSSION 

In this paper we have explored a novel idea to perform 
spectral shearing interferometry of photoelectrons using 
two coherent XUV (or x-ray) probe fields. Due to the ex¬ 
cess or shortage of photon energy for a given ionic tran¬ 
sition, the photoelectron will shift up or down in coinci¬ 
dence with the transition in the ion. The idea is closely 
related to the attosecond streak-camera method 
where a strong IR field is used to drive the electron in 
the continuum and to the PROOF method [23, 40] where 
a single IR photon is absorbed or emitted to shear the 
photoelectron distribution. The corresponding atomic re¬ 
sponse times of the attosecond streak camera are shown 
for reference in Fig. [4] (c) |38j. As can be observed, the 
response time from the outer-core method ( 2 p -> 2s tran¬ 
sition) is larger than the response time of the streak cam¬ 
era from the 2 p state in neon. As we explained, the rel¬ 
atively large response of the outer-core method comes 
from stimulated continuum electron transitions by the 
probe field. In contrast, the response time of the inner- 
core method ( 2 p -> Is) is found to be comparable to that 
of the streak-camera method. In this case the response 
time of the inner-core method comes from correlation ef¬ 
fects and possibly field-convolution effects. In theory, 
this establishes the newly proposed scheme as an all- 
XUV or x-ray method for direct group-delay determina¬ 
tion of attosecond pulses. Assuming that the attosecond 


pulse has been readily characterized, e.g. by the streak- 
camera method, the new method can be used to study 
the phase difference of two-photon (XUV-XUV or XUV- 
x-ray) processes. However, in order to extract the desired 
signal, i.e. the ^ 32 -dependent modulations of the (S+) 
peak in Fig. [2j we need to study channel-resolved photo¬ 
electrons. More precisely, we need to distinguish between 
electrons from the unexcited ion (with a 2 p hole) and the 
excited ion (with a Is or 2s hole). In practice, this is a 
major drawback of the new method because the streak- 
camera does not require any form of coincidence detec¬ 
tion. The first coincidence detection schemes combined 
with attosecond pulses have been reported recently 01 V 
143] , but so far no experiments have been reported where 
the state of the ion has been determined separately from 
the electron. Nonetheless, let us now speculate as to how 
this type of measurements can be performed in future 
experiments, inspired by the existing technology, such as 
reaction microscopes [44] and photoelectron-fluorescence 
coincidence detection [45]. First, for the case of a Is hole, 
high-energy Auger emission will occur on a femtosecond 
time scale and efficiently convert the singly charged ion 
to the doubly charged ion. A reaction microscope can be 
used to separate the photoelectron and the ion in space, 
then the ionic charge can be determined by accellerating 
the ions in an electric field. Since Auger emission is the 
dominant decay mechanism for the Is hole this Auger- 
based method is deemed more feasible than fluorescence- 
based detection. In contrast, the decay occurs exclu¬ 
sively by fluorescence for the case of a 2 s hole and it 
may appear that the only way to probe the state of the 
ion would be to detect florescence photons on a nanosec¬ 
ond time scale. However, recent experimental work has 
shown that it is possible to laser-enable Auger decay of 
the 2s hole by hitting the excited ion with an intense IR 
laser field m- This opens up for both ion acceleration 
technique and electron-coincidence detection of the high- 
energy primary electron and the low-energy Auger elec¬ 
tron to determine the state of the ion. Clearly, all these 
ideas are more challenging to implement experimentally 
than the conventional attosecond streak camera, but we 
believe that these are technical challenges that can be 
overcome in the future. Finally, we stress that the issue 
of photoelectrons with the same final energy from differ¬ 
ent states of the ion is inherent to the broad bandwidth 
of the pump pulse. If our aim is to study the phase of the 
two-photon matrix elements it is more efficient to replace 
the isolated pump pulse [(1) in Fig. [l] by an attosecond 
pulse train that translates to a comb-like photoelectron 
spectrum with spacing 25u. In this case the (S+) sig¬ 
nal then resides on peaks in between the comb-like peaks 
of the pump field, which means that the ^ 32 -dependence 
can also be studied without coincidence detection, but 
only at discrete energy intervals determined by the comb 
structure. This setup bears great resemblance with and 
could be used together with the RABBITT method [28]. 



V. CONCLUSIONS 

In conclusion, we have proposed a new type of pump- 
probe scheme that relies on stimulated core-valence tran¬ 
sitions by two narrow-band detuned XUV or x-ray probe 
fields and a short XUV pump pulse. Here, we applied 
the method to the characterization of isolated attosecond 
pulses and we demonstrated the existence of an atomic 
response time that gives insight into the nature of the 
stimulated core-valence transitions. In particular, for 
the stimulated 2 p -> 2s hole transition in neon, we found 
that the response time can be approximated by a ratio 
between electron continuum transitions and the stimu¬ 
lated hole transition. In practice, the method relies on 
coincidence detection of electron and ion, which makes it 
less efficient than existing techniques based on IR sources 
for pulse characterization. Nonetheless, the method is 
a natural candidate for future XUV-XUV experiments 
on table-top HHG sources and at FEL facilities, as it 
presents a way to study XUV-XUV/x-ray processes with 
short pump pulse probed by the sharp frequency band¬ 
width of the probe fields. 

APPENDIX A 

Time-dependent configuration interaction singles (TD- 
CIS) [35] includes the Hartree-Fock ground state |^>o) and 
its single excitations |<F^) based on the one-particle Fock 
operator H 0 and its eigenstate | ipi) with energy e^. Gen¬ 
erally, indices a,6,c... are used for spatial orbitals that 
are occupied in |<Fo), f° r unoccupied (virtual) orbitals in¬ 
dices p,g,r,...are employed, and the indices ij,/c,... are for 
general orbitals (occupied or unoccupied). Spin-orbit in¬ 
teraction is not considered in this work. The many-body 
wave packet in the CIS basis is given by 

|*,t)=ao(t)l*o) + Ei:«£(OI®S), (14) 

p a 


where z ba = {(p b \z\(p a ), P is the projection operator acting 
on the subspace composed of the virtual orbitals 

P = Z\^m = ^-Z\M(U (is) 

p a 

and J ba and K ba are, respectively, generalized Coulomb 
and Exchange operators associated with the direct 
Coulomb matrix elements v pbqa and the exchange matrix 
elements v pbaq : 

v pbqa = {(Pp\Jba\(fq) 

Vpbaq — {^plP-bal^q) - 

This procedure establishes a system of linear, coupled 
one-particle Schrodinger-like equations in Eq. ( p~7| ) for 
the orbitals \Xa(t)) with initial condition |xa(d)7/ = 0. 
To calculate photoelectron spectra, we need to evalu¬ 
ate with the transition amplitude between a modified 
Volkov state and the outgoing multi-channel wave packet, 
( Xk , a (t) | Xa(t) ), at a time, £, long after all interac¬ 
tions have ceased. The modified Volkov state in length 
gauge, |xk,a(^))? with momentum k includes an addi¬ 
tional phase-factor for a hole at the orbital |(^ a ), and it 
satisfies the equation 

*^|Xk,a(i)) = ("]V 2 -£a- E(t)z j IXk.aO) 

= H a (t) |Xk,aW), (20) 

where H a (t) represents the modified Volkov Hamiltonian 
for the electron moving in the external laser field with a 
hole fixed at the orbital \(p a )- To overcome the difficulty 
of the calculation with a large box, we adapt the time- 
dependent surface flux (t-SURFF) |36j method to the 
multi-channel TDCIS formalism. First, we define the 
overlap from a large distance R c to infinity between the 
modified Volkov state and the wave packet in a given 
channel, a, by a stepfunction: 


with initial conditions <ao(to) = 1 and o^(t o) = 0. To de¬ 
scribe the hole dynamics and the corresponding electron 
wave packet propagating in the real space, we introduce 
time-dependent orbitals that collect all single excitations 
originating from the occupied orbitals |(^ a ), 

\Xa(t)) = Y, a a( t )\ l Pp)- ( 15 ) 

P 

For the atomic system interacting with laser field E(t) 
linearly polarized along the z axis, the TDCIS equations 
of motion can be written as 

id 0 =-E(t)Y l ( l Pa\z\Xa(t)) (16) 

a 


A^ a (R c ,t) = ( Xk,o(i) | 0(Rc) | Xa(t) ) 

= [. d(3) r xi, a ( r P)Xa(r,t), (21) 

J\r\>R c 


that converges to the transition amplitude after some suf¬ 
ficiently large time T c . Eq. (21) can also be written as 


A k ,a(Ec,T c ) = [ Tc dP( XK a(tMR c )\Xa(t)) 

Jt 0 at 

= ^[(^(Xk.aC)l)^(^c)IXaO) 

+ (x k AtMR c )j t \Xa(t))\ (22) 


i^lXa(O) = (H0~£a)\Xa(t)) + £ P{K b a ~ Jba}\Xb(t)) 

- E(t)Pz{a 0 \<p a ) + |xaC))} + E (t) y z ba \xb(t)), 

b 

( 17 ) 


If we neglect correlation effects and ionic potential out¬ 
side R C1 then the time-dependent wave packet | Xa(t)) 
follows the equation of motion 

Xa(t)) = H a (t)\Xa(t)) + E(t) y Z ba \Xb(t)). (23) 
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In addition to the modified Volkov Hamiltonian, there 
is other term which makes different channels coupled to 
each other by laser field. With Eq. (20), (22) and (23), 
we can convert A^ a (R c ,T c ) from the spatial integration 
at T c to the temporal integration: 


A^ a (R c , T c ) = i £ dt[(xk,a(t)\H a (.t)6(R c )\Xa(.t)) 

- E(Xk ,a(tMRc)(H a (t)5 ab + E(t)z ab )\ X b(t))] 

= i [ dt(x\L,a{t)\[H a {t) 1 9{R c )]\xa{t)) 

Jto 

-iY, Z ab f dtE(t)(xk,a(t)\0(R c )\xb(t)) 

b Jt ° 

= - f dtJ kt a(Rc,t) - f dtK Ka (R c ,t). (24) 

Jto Jto 


with parent ions. The t-SURFF integral equation 

A^ a (Rc,T c ) = - f dtJ Ka (R c ,t) ~ f dtK Ka (R c ,t) 

Jto Jto 

(27) 

is numerically evaluated with y a (r, t) determined by TD- 
CIS under the initial condition A^ a (R c ,to) = 0. Finally, 
the momentum spectrum <7k, a (k) and energy spectrum 
&E,a(E) for the channel a are given by 

a Ka (k) = \A k , a (R c ,T c )\ 2 3 (28) 


<TE,a(E)= E 

\k\=V2E 


Ck(k) 

|k| ' 


(29) 


The numerical results presented in the main text are ob¬ 
tained using an ID effective central potential and effec¬ 
tive electron-electron repulsive potential 


We get two terms and the first term is the time integra¬ 
tion of the flux 


Jk,a( R c>t) = ^r[-xi \,a( r ^)dvXa{r,t) 

+ Xa(r,t)drXk,a( r ’ t '>] 


(25) 


Rc 


through the boundary R c from to to T c as indicated in 
E 23 EZ]. Compared with the previous work mi the sec¬ 
ond term is new and its integrand 


K Ka (R c ,t)=iY i z ab E(t)( X]i , b (t)\0(Rc)\Xb(t))e- i ^-^ t 

(26) 

represents the channel-coupling of the TDCIS via laser 
field after the electronic wave packets pass though R c . 
This contribution is missing in the integration of the flux 
at R c , so this channel-coupling term can be viewed as an 
external source from other channels as the states of the 
ion makes transition. In other words, if the field-driven 
transition between two different ionic states plays the 
role in the physical process, this channel-coupling term 
cannot be neglected. This term is especially important 
if the photoelectron spectrum is measured in coincidence 


Ves(z) 


ZeS 

\A 2 + 2c 


kee (% 11 % 2 ) — 


y(^i-z 2 ) 2 +^ 


(30) 


with parameters Z e ff,Z ee ,z c , and z e that reproduce the 
experimental ionization energies of the 2s and 2p or¬ 
bitals in neon. We considered two models: First, the 
parameters of the effective potential were chosen such 
that the electron-electron interaction was zero, which 
corresponds to the independent particle approximation 
(IPA), by parametrization as Z ee = 0, Z g q = 1.795, and 
z c = 0.7. Second, we studied the correlated TDCIS 
model, parametrized by Z ee = 1, Z e ^ = 9, and z c = 0.755. 
The IPA and TDCIS agree remarkably well for our field 
parameters with only slightly different extracted phase 
parameters for the proposed method that depend on the 
detailed correlated interactions. 
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